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^ ■ Abstract 

' Event generators are an indispensable tool for the preparation and analysis of 

particle-physics experiments. In this contribution, physics principles underlying the 
^ , construction of such computer programs are discussed. Results, within and beyond 

I the Standard Model of particle physics, obtained with a new event generator [1] are 

■ presented. This generator is capable to describe signal processes for exotic physics 

, and their backgrounds at electron-positron and proton-(anti)proton colliders. 
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^ ■ 1 Introduction 



In current collider experiments the hard interaction of two incoming beams 
results in the production of up to thousands of outgoing particles. So far, no 
evidence has been found that contradicts the belief that this process is de- 
scribed by the Standard Model (SM) for strong and electroweak phenomena. 
Unfortunately, a full quantum-mechanical treatment is out of reach. There are 
two reasons for this apparent shortcoming: first of all, the sheer number of par- 
ticles involved gives rise to a tremendous number of interfering contributions 
that grows factorially with the number of particles. Furthermore, perturbation 
theory is not able to account for the transition of partons to hadrons. This 
failure of perturbation theory necessitates other methods. Simulation tools 
called Monte Carlo programs or event generators have proved to be successful 
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Fig. 1. Sketch of a proton-proton collision from an event-generator point of view. 

for a detailed description of multiparticle production. In this contribution we 
want to briefly illuminate the principles underlying such programs, and we 
would like to present some results of our new event generator SHERPA [1]. 



2 Physics in Event Generators 



For a brief discussion of the physics concepts underlying event generators, 
Fig. 1 will be employed, in which a proton-(anti)proton collision is sketched. 
The simulation of such an event has to be decomposed into different indepen- 
dent stages. To control this independence, these stages are ordered according 
to the scales related to them. 

At the highest scale the signal process takes place leading eventually to the 
production of heavy and/or exotic states. In the figure, this is represented 
by the thick blob in the middle, in which two gluons fuse to produce a four- 
particle final state. For such signals one would like to use the full matrix ele- 
ment. Within our event generator the module AMEGIC++ [2] takes care of the 
hard processes. It automatically creates all Feynman diagrams, stores them as 
helicity amplitudes into library files and uses these libraries when the integral 
over the phase space of final-state particles, i.e. the cross section, is evaluated 
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through Monte Carlo methods. 

As the event evolves down to lower energy scales, multiple emission of sec- 
ondary particles must be described. The bulk of radiation activity is in the 
soft and/or collinear phase-space region, where the interference of subsequent 
emissions is of subleading logarithmic order only. The logarithms are in ra- 
tios of the actual parton-emission scale compared to a suitably defined harder 
scale, so that, for a description of the radiation pattern, these interferences can 
be safely neglected. Hence multiple parton emission factorizes into chains of 
individual parton branchings. Consequently, one obtains a simulation in terms 
of a parton shower. Radiation off the initial state is more involved, but must 
result in a boost only, which, due to recoil effects, may also produce some net 
transverse momentum of the final state. In our event generator, the parton 
showers are performed by the module APACIC++ [3], which resums all leading 
logarithms and uses the invariant masses of the decaying partons to order 
subsequent emissions. An ordering in successive opening angles to properly 
account for quantum coherence effects is supplemented. A novel feature of our 
code is that, for the first time, an algorithm has been developed combining 
matrix elements and the parton shower by smoothly filling the parton-emission 
phase space and by respecting the logarithmic accuracy [4]. 
Some particles produced in the hard interaction might be unstable. If their 
mass is large enough, their decay might involve particles emitted with large 
energies or under large opening angles. Such decays, again, are treated by 
using matrix elements and their subsequent parton showers. At high energies 
the possibly perturbative interaction of the beam remnants, sometimes also 
coined multiple parton interaction, comes into play. These interactions may 
develop their own parton showers. However, this phenomenon has the poten- 
tial to severely obscure the interpretation of experimental data and, therefore, 
has to be simulated. These two aspects are currently being implemented. 
The gauge structure of QCD leads to a breakdown of perturbation theory at 
scales of the order of Aqcd- So, the parton shower has to stop at such scales 
and its parton ensemble has to be translated into hadrons, i.e. white objects. 
Since confinement has not been completely understood yet, hadronization ef- 
fects have to be modelled phenomenologically introducing parameters that 
have to be fitted to data. Successful models, such as the string and the cluster 
model, rely on the fact that before the non-perturbative transition the parton 
shower has a definite colour structure ^ . Our code incorporates a newly de- 
veloped cluster-hadronization model and an interface to the well-known Lund 
string fragmentation, allowing for systematic studies of model uncertainties. 



The parton shower is inherently in a large- A'^c approximation, where A'^c is the 
number of colours. Employing a forced splitting of gluons into qq pairs, the partons 
form a set of colour singlets, each of which contains a quark and an antiquark. 
In cluster models usually formulated in terms of formation and decay of clusters 
these singlets are interpreted as primordial hadron resonances with indefinite mass 
decaying into primary, specified hadrons. 
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3 Results 



3.1 Results obtained with AMEGIC++: the parton level 

Here, we will highlight the abilities of AMEGIC++ in calculating cross sections 
for processes within and beyond the SM. There is a plethora of reasons to 
search for new physics, e.g.: 

The SM has two intrinsic problems, the hierarchy and the naturalness problem, 
which are addressed to the question why the scale of electroweak symmetry 
breaking and the Planck scale, where strong gravitational effects are expected, 
are separated by a wide desert. To understand this hierarchy self-interactions 
of the (scalar) Higgs boson are considered. Using a cut-off procedure for the 
ultraviolet divergencies with the cut-off given by the highest scale in the SM, 
the Planck mass, one immediately realizes that the counter term for the Higgs 
mass has to be quadratic in the Planck mass. In reality, however, unitarity of 
WW scattering forces the Higgs mass to be lower than roughly 1 TeV. The 
difference of roughly seventeen orders of magnitude requiring an extreme fine- 
tuning (naturalness problem) lets the SM seem somewhat inconsistent. 
On top of this intrinsic problems, gravitation is not incorporated by the SM. 
Assuming the existence of a "theory of everything" , one may conclude that 
the SM is incomplete and constitutes only an "effective theory" . 
However, the hope to successfully search for new physics depends on our ca- 
pabilities to understand and subtract the "old" physics, i.e. the SM. Thus it is 
crucial for any event generator to provide an optimal simulation of SM effects. 

3.1.1 Multi-particle production in the SM: six-fermion final states 

Despite that six-fermion final states constitute significant backgrounds to sig- 
natures for new physics, they provide unique opportunities to study some of 
the most interesting aspects of the SM in great detail. The most important 
channels are the production and decay of pairs of top quarks and - if existent 
- of one or more Higgs bosons, the latter allowing to test the structure of the 
Higgs potential. If no evidence for a Higgs boson was found at the LHC, the 
study of quartic gauge-boson couplings is mandatory in order to understand 
alternative scenarios of electroweak symmetry breaking. To test these complex 
calculations involving up to (9(10^) Feynman diagrams, tuned comparisons of 
different generators are required. Table 1 presents predictions of LUSIFER [5], 
HELAC/PHEGAS [6] and AMEGIC++ [2] for some six-fermion final states corre- 
sponding to the production and decay of a top-quark pair. More extensive 
comparisons of results obtained by HELAC/PHEGAS and AMEGIC++ also taking 
into account finite fermion masses can be found in [7]. 

3.1.2 Cross sections for Supersymmetry 

The ability of our code for searches of new physics is exemplified by some 
processes within the framework of the MSSM - the Minimal Supersymmet- 
ric extension of the Standard Model, see e.g. [8]. Supersymmetry is one of 
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a[fb] 


LUSIFER 


HELAC/PHEGAS 


AMEGIC++ 




5.853(7) 


5.866(9) 


5.879(8) 




5.819(5) 


5.822(7) 


5.827(4) 




5.809(5) 


5.809(5) 


5.809(5) 


u^^'^T~ i>T-hh 


5.800(4) 


5.798(4) 


5.800(3) 




17.171(24) 


17.204(18) 


17.209(9) 



Comparison of tree-level predictions corresponding to six-fermion top-quark pair- 
production channels at Ecm = 0.5 TeV for massless fermions; see [9]. 

the most promising candidates for physics beyond the SM, since it cures the 
hierarchy problem in a natural fashion: it predicts at least one superpartner 
for each particle, a bosonic (fermionic) partner for a fermion (boson), with 
identical quantum numbers and mass. The only difference is in the spins, 
which differ by one half. In loops, the partners come with opposite sign, thus 
the quadratic divergencies mentioned before are cancelled. However, no su- 
perpartner has been found yet. Hence, Supersymmetry must be broken. The 
demand to solve the hierarchy problem implies a soft breaking, i.e. the form 
of breaking terms is limited. First, one usually considers a minimal Supersym- 
metry, i.e. one superpartner per particle of the SM, and this MSSM is broken 
by adding all possible breaking terms. 

The complete MSSM has been implemented in AMEGIC++, and an example is 
presented in Fig. 2: gaugino-pair production at an e~^e~ collider shows a strong 
dependence on the masses of the sleptons propagating in the t-channel. This 
is true for the lightest charginos (xt) and neutralinos (x?)- Starting with the 
mSUGRA point ^ SPSla [10], the appropriate slepton masses have been var- 
ied. For the charginos, the rise for large sneutrino mass exhibits the vanishing 
destructive interference of s- and t-channel diagrams. 



Neutralino pair production at E^^^^ = 500 GeV Chargino pair production at E^^^^j, = 500 GeV 

0.4 0.5 I — • — • — • — ■ — , — • — • — • — ■ — , — • — • — • — • — , — • — • — r- 




qI 0.1 1 

100 150 200 250 300 350 400 100 150 200 250 300 

physical selectron mass (GeV) physical sneutrino mass (GeV) 

Fig. 2. Ecm = 500 GeV cross sections for e~^e~ — > XiXi (left) and for e+e" —>■ xfxi 
(right) in dependence on the selectron and the sneutrino mass, respectively. 

^ For current studies of the phenomenology of the MSSM, specific benchmark points 
have been proposed [10]. 
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Besides e~^e~ physics, a next linear collider might provide an unique opportu- 
nity to study 77 and 76 interactions at high luminosities. Both can be sim- 
ulated within our event generator. High-energetic electrons are transformed 
through Compton backscattering of laser light into high-energy photons. This 
energy spectrum is shown in the left part of Fig. 3. Using the polarization 
dependence of the backscattered photons, valuable information on particles 
produced in collisions of polarized beams can be obtained. The right part of 
Fig. 3 depicts the photon polarization in dependence on its energy-fraction 
w.r.t. the incoming electron. Figure 4 illuminates smuon pair production ex- 
hibiting a significant dependence on the polarization states of the 7 beams. 




Fig. 3. The photon luminosity spectrum and the degree of circular polarization as 
a function of the photon energy-fraction x = Ey/Eg for different laser and beam 
helicities in accordance to the parametrization by CompAZ [11], assuming an electron 
energy E(, of 250 GeV and a laser energy u>j of 1.17 • 10^^ GeV. 



0,2 



0,15 



O 0,1- 



0,05 





'1 














\-+ 








H- 

00 
















,1 




= 100 GeV 



200 



400 



600 
Ecm/G'=V 



800 



1000 



Fig. 4. The total cross section of 77 — fl'^fi for L = (++), L = 2 (+- 
completely unpolarized beams (00) using the CompAZ spectrum, note S = /i. 
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3.1.3 Cross sections for Large Extra Dimensions 

An alternative physics scenario has been recently put forward by Arkani- 
Hamed, Dimopoulos, and Dvali (ADD model) [12]. They propose to solve the 
hierarchy problem of the Standard Model by shifting the Planck scale down to 
scales of the order of TeV. This is achieved by introducing n > 2 new compact 
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Fig. 5. a) Differential cross section for e~^e~ . Tlie solid- line histogram 

shows the SM background. Crosses and circles depict virtual graviton exchange in 
two ADD scenarios, b) Total cross section for e'^e~ jG as a function of -Ecm- 

spatial dimensions that are large compared to the scale of electroweak sym- 
metry breaking. These dimensions can in principle be as large as fractions of 
millimeters, since only down to these length scales gravity has been probed ex- 
perimentally. To ensure that the SM remains valid within the regions tested so 
far, SM particles remain constrained to the common four-dimensional man- 
ifold. Only gravitons propagate in all dimension. This changes the laws of 
gravity and leads to a number of striking signals for the next generation of 
collider experiments. To study the rich phenomenology of this ADD model, 
a correct treatment of spin-2 particles in the helicity formalism has been de- 
veloped and implemented [13]. This allows to describe real as well as virtual 
graviton production, both of which exemplified in Fig. 5. 

3.2 Results obtained at the hadron level 

They are presented for e~^e~ collisions at the pole, Ecm = 91-2 GeV. We 
will focus on the overall performance of our event generator in its ability to 
describe multi-hadron production, and on some aspects of our cluster model. 

3. 2. 1 Overall performance 

First, inclusive observables are considered. This is usually done using event- 
shape variables, such as thrust, sphericity, oblateness etc., or by looking into 
the correlations of hadronic matter - jets - deposed in the detector, such as 
relative angles of jets. In Fig. 6 oblateness and 1— thrust are depicted. The 
results of our new generator SHERPA [1] are confronted with the results of the 
well-known PYTHIA Monte-Carlo [14] and with experimental data obtained 
by ALEPH [15]. Both predictions agree well with the data. Moreover, using 
DELPHI data [16], Fig. 7 shows the comparison for two more exclusive ob- 
servables, the and the Bengtsson-Zerwas relative angle in four-jet events. 
There, the impact of our new approach to merge matrix elements and parton 
showers proves its predictive power. In contrast, PYTHIA adds radiation to the 
initial two-jet configuration mainly through the parton shower. Therefore, it 
includes less quantum interference effects. This is the source of the deviation. 
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Fig. 6. Event-shape variables for e~^e — > hadrons at the pole. 
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3.2.2 Cluster-hadronization model 

Inspired by the well established Webber model [18] for primary hadronization 
of a coloured partonic system, we have developed a new model presented in 
[19] and extended by the following features: soft colour reconnection is mod- 
elled and incorporated in the formation and fragmentation of clusters, yielding 
additional cluster-decay configurations. A new mesonic cluster type, the two- 
diquark cluster, and its decay channels are added. The explicit consideration 
of the spin of diquark cluster-constituents is rendered possible. The most sig- 
nificant feature of our treatment is the flavour-dependent separation of the 
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Charged particle log( 1 / x distribution 
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Fig. 8. = ln(l/xp ) distribution of charged particles for Ecm = 91.2 GeV 
in e'*'e~ annihilation into light quarks only. The preliminary SHERPA predic- 
tion is presented together with OPAL data [20], and with results from default 
PYTHIA-6 . 1 {uds) and default HERWIG-6 . 1 {uds) [21]. 

cluster and hadron regimes in terms of the cluster's mass. As a consequence, 
the different cluster-transition modes are obtained automatically. 
A preliminary SHERPA version including our cluster-hadronization model is 
capable to describe e^e~ ^* jZ^ — * dd,uu,ss —>■ hadrons. The agreement 
with default PYTHIA-6 . Knds) [14] and experimental light-quark data is sat- 
isfactory. To illustrate this, the distribution of the negative logarithm of the 
scaled-momentum, ^^'^'^ = — logx^'^'^, is shown in Fig. 8. 



4 Conclusion 

Event generators will continue to be an indispensable tool for future exper- 
imental analyses. A new multi-purpose generator for high-energy collisions, 
SHERPA [1], is currently being developed. It already proved to successfully de- 
scribe e+e~ annihilation. The embedded automated matrix-element calculator 
AMEGIC++ [2] is not only capable to evaluate all Standard Model processes, but 
also physics beyond the Standard Model (e.g. MSSM and ADD). 
The code is written in a transparent structure using the object-oriented pro- 
gramming language C++. Consequently, the implementation of extensions is 
straightforward; the code is easy to maintain and adjust to the needs of users. 

We acknowledge that this work has been supported by BMBF and GSI. 
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